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NAVAL ARCHITECTURE 


Manoeuvring forces on azimuthing podded 
propulsor model 


Maciej Reichel 
Ship Design and Research Centre 
(CTO S.A.) Gdansk 


ABSTRACT 


This paper presents the preliminary part of comprehensive manoeuvring open- water tests 
of a gas carrier model. The paper focuses on open water experiments with an azimuthing 
podded propulsor. The test program was carried out in the cavitation tunnel and the large 
towing tank of Ship Hydromechanics Division, Ship Design and Research Centre, Gdansk. 
The pod was tested as a pushing unit with a 161.3 mm diameter propeller. Steering forces 
were measured in the range of advance coefficient from 0.0 to 0.8 combined with the range 
of deflection angles from -45° up to +45°. Measurements on the pod without propeller 
were also performed. The experiment results are presented in the form of non-dimensional 


coefficients in function of advance coefficient and deflection angle. Analysis of the experimental results 
and the conclusions are presented. 


Keywords : podded propulsors, open water experiments, manoeuvring forces. 


INTRODUCTION 


Podded propulsion, an integrated propulsion and steering 
system is known for nearly half a century [1].This type of pro- 
pulsion has many advantages in comparison with traditional 
propeller-rudder set, e.g. vibration reduction and flexibility of 
internal design arrangement. But the most important advantage 
is its better manoeuvrability. 

Although the good manoeuvring abilities of podded ships 
are known the manoeuvring characteristics are still not fully 
recognized. Standard experiments with free running propel- 
lers are not sufficient because the interaction effects between 
propeller and pod housing are very difficult to predict. The- 
refore it is advisable to carry out model tests with different 
podded drives to get data basis for the forces and moments at 
different steering angles [2, 3]. It is important to collect data 
for many possible pod configurations, e.g. push, pull, single 
or twin set. 

The knowledge of the forces and moments acting on the 
propeller and podded drive at different steering angles is ne- 
cessary for design and optimisation of azimuthing propulsion 
systems. Also information about wake velocity field is impor- 
tant in design aft part of the hull and possibly also skegs or 
additional rudders [4]. 


MODEL OF THE PODDED DRIVE 


Tab. 1. Main data of P447 propeller model . 


Propeller diameter D mm 161.3 
Pitch at r/R = 0.7 P mm 126.0 
Pitch ratio at r/R = 0.7 P/D - 0.785 


Expanded blade area AJA i 0.550 
ratio K 
Hub ratio d,/D - 0.219 
Number of blades Z - 4 


Direction of rotation 


Revolutions per minute 


Gondola length 
Gondola diameter 


Strut height (from 
propeller shaft) 


Strut length 


Fig. 1. Arrangement of pod propulsors and pod geometry . 
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TEST FACILITY AND PLAN 
OF THE EXPERIMENT 


The open water tests were carried out in the cavitation tunnel 
and in the large towing tank of Ship Hydromechanics Division, 
Ship Design and Research Centre in Gdansk, (OHO-CTO). 


Fig. 2. Podded drive in the OHO-CTO large towing tank before tests . 


Measurements were conducted at the constant number of 
revolutions equal to 900 rpm. The carriage velocities were so 
established as to obtain advance coefficient of 0.2; 0.4; 0.6; 0.8. 
The experiments were performed for thirteen deflection angles 
in the range from -45°up to +45°, with the step of 5° from -15°to 
+15°and 10 within the rest of the range. Forces on the bare pod, 
i.e. without propeller, were also measured. The positive deflec- 

tion angle and directions of forces are shown in Fig.3. 


Fig. 3. Definition of forces . 


RESULTS 


All the results are presented in the traditional 
non-dimensional form by applying the following formulae : 


U 
= Advance coefficient J= a (1) 
= Longitudinal force coefficient FX = pD (2) 
Fy 
= Transverse force coefficient Ky z DnD? (3) 
pn 
| | 
= Thrust coefficient T pn2D4 (4) 
; _ N 
= Normal force coefficient K= a4. (5) 
pn‘ D 
| k- Q 
= Torque coefficient Q- pnD4 (6) 
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where : 


U „s inflow velocity 

ô -deflection angle 

p -water density 

n -propeller revolutions 
D -propeller diameter. 


The thrust and normal force were calculated by applying the 
following transformation from the towing tank coordinate 
system to the propulsor coordinate system : 


T|_|cosd -sinô | Fx 
N| |sind cosd | Fy 
In Fig.4 the open water characteristics for the P447 propeller 


are shown. The characteristics were measured on the model 
without pod housing and only for straight flow. 


0.70 
No - d 
0.60 7G S 
10 KQ o- 7 


KT 7 x 


0.50 


0.40 7 


0.30 


0.20 


0.10 


0.004 
0.00 0.10 0.20 0.30 040 0.50 0.60 0.70 
Fig. 4. Open water characteristics of P447 propeller . 


Results of the measurements of the forces on the bare pod 
are presented in the figures below. The experiment was made 
for one advance coefficient equal to 0.6 and for five deflection 
angle values : +45°, +10°, 0°. 

- 0.00- 


- -0,07- 


0.04- 


50 


l i O04: : l 
Fig. 6. Transverse force coefficients vs. deflection angle of bare pod . 


- 0.08- 


Fig. 8. Normal force coefficients vs. deflection angle of bare pod . 


The subsequent figures show results for the podded 
drive with right-handed P447 propeller. 


0.40 


0.35 
0.30} 
0.25} 
0.20} 
0.15 
0.10 
0.05 l- 


0.00 


-0.05 


pod angle [deg] 


-0.10 
-50 -30 -10 0 10 20 30 40 50 
—a— 0.00  —e 0.20 -—*- 040 #060 —*—0.79 


Fig. 9. Longitudinal force coefficient vs. deflection angle . 


* the longitudinal force coefficients show parabolic dependen- 
ce on deflection angle, with the tendency to reduction for 
higher angles. For positive deflection angles the reduction 
of longitudinal force coefficients is higher than for negative 

ones (Fig.9) 


* the transverse force coefficients are approximately linear 
in function of deflection angles, with deviation to higher 
values for larger deflection angles and advance coefficients. 
In the vicinity of 15° deflection angle a little disturbance is 

observed (Fig. 10) 


* the thrust coefficients depend strongly on deflection angles 
and propeller loading. For both steering directions the thrust 
is increasing but its increase is visibly stronger for negative 

deflection angles (Fig.11) 


* the normal force coefficients behave like the transverse 
force coefficients, i.e. approximately linear in function of 


0.20 


-35 -30 -25 -20 -15 -10 -5 0 
Pod angle [deg] 
—s— 0.40 


—o 0.60 —— 0.79 


20 25 30 35 40 
Pod angle [deg] 


Fig. 10. Transverse force coefficient vs. deflection angle . 


0 5 10 15 


-20 -10 0 10 20 30 40 50 
Pod angle [deg] 


—a— 0.40 


—e 0.20 —e— 0.60 —— 0.79 


Fig. 11. Thrust coefficient vs. deflection angle . 
deflection angle, however with a faster increase for higher 
advance coefficients. Like in the case of the transverse force 


coefficients a little disturbance is visible in the vicinity of 
15° deflection angle (Fig. 12). 


To calculate the coefficients of forces in function 
of advance coefficient the following formulae were applied : 
Uer 7 Ucar 6080 (7) 
where : 
U 


* for smaller pod deflection angles the longitudinal force 
coefficients are decreasing faster than for higher ones. 
The nature of the dependence is similar to that of a classic 

propeller (Fig.13) 


car ~ towing carriage velocity. 


* for small pod deflection angles up to 15°, the transverse 
force coefficients depend almost linearly on advance co- 
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0.30 


0.25 


0.20 


0.15 


-35 -30 -25 -20 -15 -10 -5 0 
Pod angle [deg] 


—o 0.20 —+— 0.40 —e— 0.60 —x— 0.79 


0 S5 10 15 20 25 30 35 40 45 50 
Pod angle [deg] 


Fig. 12. Normal force coefficient vs. deflection angle . 


$ $ $ $ : JE] 
0.20 0.30 040 0.50 0.60 0.70 0.80 0.90 


0.20 030 040 0.50 0.60 0.70 0.80 0.90 


Fig. 13. Longitudinal force coefficient vs. advance coefficient . 


efficient; for higher angles the dependence is changing to 
parabolic one (Fig. 14) 


* the thrust coefficients in function of advance coefficient are 
decreasing but for positive pod deflection angles the decre- 
asing is independent on deflection angle value (Fig.15) 
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JH 


-0.05 i x : j j z Pi 

0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 
0.05 

KFy [-] 


-0.30 
0.00 0.10 


0.05 


0.00 i i 
0.00 0.10 0.2 


0.05 oes eee ees eenees: : sesesess 


0.00 i i 
0.00 0.10 0.60 


Fig. 15. Thrust coefficient vs. advance coefficient . 


0.70 0.80 0.90 


0.20 0.30 0.40 0.50 


* the normal force coefficients show tendency to change their 
character from linear to parabolic for higher pod deflection 
angles in function of advance coefficient (Fig.16) 


Fig.17 introduced to summarize the content of the previo- 
us figures shows proportion of thrust force and normal force 


between positive and negative pod angle. The value of 100% 

represents the thrust force and normal force for the positive 

direction of pod angle. For the negative pod deflection angles 
the absolute value of normal force was used. 


TE 


0.60 0.70 0.80 0.90 


0.60 


Fig. 16. Normal force coefficient vs. advance coefficient . 


0.20 0.30 0.40 0.50 0.70 0.80 0.90 


KNdelta+/ 


0 5 10 15 20 25 30 35 40 45 50 
Deflection angle [deg] 
—+— 0.40 —e— 0.60 —— 0.79 


—a— 0.00 —e— 0.20 


160 


: KTdėlta+/ : 


140 
120 
100 


0 5 10 15 20 25 30 35 40 45 50 
Deflection angle [deg] 


Fig. 17. Interdependence between KT and KN for positive 
and negative deflection angles 


* the ratio of normal force coefficients shows difficult depend- 
ence for small deflection angles, for angles higher than 15° 


the relationship shows tendency to be linear in the vicinity 


of in the range of 80 — 100 % AS we . The exception to the 
Se 

tendencies is visible at J = 0.00. For small advance coef- 

ficients the normal force coefficient ratio is smaller than 
for higher ones 


* the ratio of the thrust coefficients shows tendency to increase 
for small pod angles up to 10°, and to decrease for higher 
deflection angles. This tendency is more visible for higher 

advance coefficients. 


CONCLUSIONS 


To summarize results obtained from the towing tank tests 
the following conclusions may be presented : 


O the presented results clearly show the typical hydrodynamic 
characteristics of azimuthing podded propulsion 


O the asymmetries in values of the force coefficients for posi- 
tive and negative deflection angles are due to the influence 
of the direction of propeller rotation (to the right) 


O with positive pod angle a negative normal force is produced 
and vice versa, which results in a destabilizing moment 
tending to increase the turn rate 


O the increase of thrust coefficients for higher deflection an- 
gles is induced by reduction of axial inflow velocity which 
involves reduction of advance coefficient 


O for deflection angles larger than 15° the thrust and normal 
force coefficients are from 10% to 30% smaller than for the 
corresponding positive deflection angles, which is caused 
due to the interaction between right — handed propeller and 

podded drive. 
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NOMENCLATURE 


A,/A, - expanded blade area ratio 
c - strut length 

d - gondola diameter 

- hub ratio 

D - propeller diameter 

Fx - longitudinal force 

Fy  - transverse force 

h - strut height (from propeller shaft) 
J - advance coefficient 

K - longitudinal force coefficient 
K - transverse force coefficient 
K - normal force coefficient 

K - torque coefficient 

K - thrust coefficient 

i - gondola length 

- revolutions per minute 

- normal force 

pitch at r/R=0.7 

pitch ratio at r/R=0.7 

- torque 

- thrust force 

carriage velocity 

- inflow velocity 

- number of blades 


wuzs 
J 


a 
S, 


- pod angle 
- water density 


Da NACHO 
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«~~ C onference — 


on Exploitation Foundations 


Maritime University, was held a scientific seminar of the 
Regional Group of the Section on Exploitation Founda- 
tions, Machine Building Committee, Polish Academy 
** On failures of piston-cylinder system of ship internal 


** Investigations of ageing process of lubricating oils in 


“* Anew concept of building ship power plant simulators 


REGIONAL GROUP 
of the Section 


On 18 January 2007 at Mechanical Faculty, Gdynia 


of Sciences (PAS). 


The Seminar program contained presentation of six pa- 
pers elaborated by scientific workers of the Faculty, 
as follows : 


œ Professor Jan Kazimierz Włodarski — 55 years of 
activity in the field : construction and operation of 
machines — by R. Cwilewicz and K. Witkowski 


combustion engines — by J. K. Wtodarski 

œ Assessment of influence of indicator channel and its 
dimensions on results of pressure measurements in 
engine cylinder — by W. Gałecki 


trunk-piston engines — by A. Mtynarczak 

œ Possibilities of diagnosing selected failures of injection 
devices of self-ignition engines on the basis of indicator 
diagram run — by R. Pawletko 


with the use of 3D visualization — by L. Tomczak 


After interesting discussion the Seminar participants 


a visited scientific laboratories of the University. J) 
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Miscellanea 


Jubilee of 50" Anniversary 


On 30 June — 1 July 2006 the Experimental Centre 
for Ship Model Tests in Open Waters, Faculty of Ocean 
Engineering and Ship Technology, Gdansk University 
of Technology, celebrated 50" anniversary of its rich 

activity. 


The Centre has been established due to initiative of 
Prof. Lech Kobylinski, Head of the then Ship Theory 
Department, Faculty of Ocean Engineering and Ship Tech- 
nology, Gdansk University of Technology. All the matter 
have started from hydrofoils. In 1955 a small, 4.5 m long, 
hydrofoil was built in the Department, whose first trials 
was conducted on Motława river, and simultaneously the 
Department was assigned to conduct realization of a pro- 
ject for development of hydrofoils, in which Polish Navy 
was interested. In that time Prof. Kobylinski decided to 
organize the Experimental Centre on the area of a deserted 
water sport centre by Jeziorak lake in Ifawa. In summer 
1956 just in this place the research trials on hydrofoils star- 
ted again. Results of the research on a few manned models 
as well as relevant theoretical elaborations resulted in the 

mastering of design principles of such floating units. 


In 1965 in Wista Shipyard was built a passenger hy- 
drofoil for 76 persons, which was operated by Szczecin 
Coastal Shipping Co till 1968, and after modernization - by 
Polish Navy on the Gdynia - Hel route till 1993. Finally the 
activity in that area was ended in the 1960s when several 
hydrofoils were purchased in Soviet Union. Apart from 
hydrofoils, model tests started in the range of resistance, 
propulsion, propellers, manoeuvrability, sea-keeping 
qualities of conventional ships, and also of river pushers 
and push-trains, hovercrafts and fast boats; also various 

special investigations were carried out there. 


Results of the tests were implemented by shipyards’ 
ship design offices, shipping companies and classifica- 
tion institutions, as well as they were presented at many 

scientific conferences in Poland and abroad. 


25 years ago, also in Iława, has started the first ship 
manoeuvring training course for ship masters, carried out 
with the use of large manned ship models, that contributed 

to arranging the Training Centre by the Silm lake. 


These services and research carried out in Hawa 
have initiated organization of symposia on ship hydro- 
mechanics in Poland, and later - also symposia on ship 

manoeuvrability. 


Hence it is not an exaggeration to state that the Centre 
in question deserves the name of the “cradle” of Polish 
ship hydromechanics. 


In present, the Centre fulfils also the role of the base 
for international student scientific camps, and permanently, 
the student training base for yachting sports. 
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On the fatigue-critical amplitude 
of random-amplitude stress 


Janusz Kolenda 
Gdansk University of Technology 
Polish Naval University 


ABSTRACT 


Uniaxial non-zero mean stress of constant circular frequency in the high-cycle fatigue 
regime is considered. It is assumed that equation of the S-N curve and modified Soder- 
berg equation are applicable. For constant-amplitude stress, the fatigue-critical stress 
amplitude is defined as that which leads to failure during the required design life. For 
random-amplitude stress, expected values of the fatigue-critical stress amplitude and total 
fatigue damage accumulated during the required design life are estimated. It is found that 
the probability of fatigue failure is equal to the probability of exceedance of the fatigue- 
-critical stress amplitude. As an example, for stationary random stress the equivalent 


random-amplitude stress and probability of fatigue failure are determined. 


Key words : uniaxial load, random stress, high-cycle fatigue, failure probability 


INTRODUCTION 
If a metallic element is subjected 
to asymmetric tensile-compressive stress 
O(t)=o,) +0, sinat (1) 


in the high-cycle fatigue regime the following equations 
are frequently applied [1, 2] 


No™=K (2) 
o, = ol- s Zu osL (3) 
where : 


K — fatigue strength coefficient 

L -— maximum stress amplitude satisfying Eq. (2) 
(above which low-cycle fatigue may occur) 

m — fatigue strength exponent 

N — number of cycles to cause fatigue failure 

R, — tensile yield strength 

Z_,— fatigue limit under fully reversed tension-com- 
pression 

o — amplitude of the fully reversed stress at a given 
number, N, of cycles to cause failure 

o, — amplitude of the stress (1) which will give that 
fatigue life 

6, — mean stress 

œ — circular frequency. 


If N, denotes the required cycle number 
the design criterion N} < N becomes : 


m 


K Oy 


R. 


I= 


for : 


Baai 


e e 


Zeler 


Hence the stress amplitude o,, that leads 
to fatigue failure in N, cycles is : 


26, | K” (5) 
R. Na 

This quantity will be called 
the fatigue-critical stress amplitude. 


DESIGN CRITERION AT STATIONARY 
RANDOM-AMPLITUDE STRESS 


O5 


Let us determine the design criterion corresponding 
to Eq. (4) in the case of uniaxial Gaussian stress : 


0, (t)= 6) +6, sin(mt + a) (6) 


with the random amplitude o, and random phase angle a. 
As to the stress amplitude it follows Rayleigh distribution [3] : 


- 2) (7) 


where : s_is the standard deviation of the stress amplitude. 


F(o,) = 1-—exp 


In terms of ensemble averages, 
the counterpart of the criterion (4) reads : 


(8) 
where : E{"} denotes the expected value. 
In the considered case : 


E{o,}=(0.51)"’s, (9) 
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Blor}=2"7r (14+ >)s? (10) 


where : T is the gamma function. 
So, the criterion (8) becomes : 


m/2 
PN, 


K(I- Co) 
for : 


Oo 1/2 
= RL) < O.S) LaL 


r(i+S)sv<i (11) 


Oo 
Z| T 
FATIGUE-CRITICAL STRESS AMPLITUDE 
AND PROBABILITY OF FATIGUE 


FAILURE UNDER STATIONARY 
RANDOM-AMPLITUDE STRESS 


Eqs (9) and (11) yield : 
NE (1+ FENH” 


„oj K 


m/2 
(0.257) ( R 


Hence the expected value of the fatigue-critical 
amplitude of the stress (6) can be estimated as : 


o) L j 13 
R, Nar (1+) i 


With the respect to the distribution function (7), the proba- 
bility, P „, ofexceedance of the fatigue-critical stress amplitude 
is given by : 


<1 (12) 


= (0.251 (1 


2 


Cor 
P,, = P{0,> 64} = exp (- 38) (14) 
From Eqs (13) and (14) one obtains : 
2 2/m 
_ 0.257 _ So K 
Po Ta exp 252 ( R ) N r ( m) (15) 
X e d + 2 


Now, the task is to determine the relation between P.,, and 
the probability, P ,, of fatigue failure under the stress (6). Accor- 
ding to the Palmgren-Miner rule [1, 2] the increment of fatigue 
damage, caused by N, stress cycles of constant amplitude is 
equal to the ratio N/N where N is the cycle number to cause 
failure. Referring to the stress (6) the expected cycle number 
to failure can be estimated by means of Eq. (12) as : 


(0.257)™?(1- So) K 


N= = (16) 
m m 
r(l+™)(Ete 3] 
Hence the expected value, D, 
of the fatigue damage after N, cycles is : 
Nr(i+™ Ee i" 
D= a (17) 
m/2 0 
(0.25n)™(1- se) K 
and the probability of fatigue failure : 
P.=P{D>1} (18) 
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becomes : 


1/m 
K 
P,=P4o,2(0.25n)'(1- So) —| |. 
e Nar (i+ 2) 
=P[0 20a] (19) 


Hereby the following has been proved : In the high-cycle 
regime the probability of fatigue failure under random- 
-amplitude stress is equal to the probability of exceedance 

of the fatigue-critical stress amplitude. 


EXAMPLE 


Let us consider a stationary random stress : 
5, (t)=o9+ 6, (t) (20) 


where : o, is the mean value of the stress 6, (t) and G,(t) is 
the stationary (in the wide sense) stochastic process of zero 
mean value and known power Spectral density S(@). The 
suggestion is to model the stress © (t) by the equivalent stress 
o (t) in the form of a periodic (in the mean- square sense [4]) 
Gaussian process : 


6,(t)=o, +asin(@,t + @) = Ql) 
=O, + a,exp(j@,t)+a,exp(— jo,t) 
where : 


6, — the mean value of the stress S(t) 

©, — the circular frequency of the equivalent stress 
a — the random amplitude of Rayleigh distribution 
ọ — the random phase and : 


a i $ 
a1 = exp UO) , a2 =a] 


E{a,} =E{a,} =E{aa;}=E{aja,}=0 


(C) — complex conjugate 
j — imaginary unity. 


Following the approach based on the theory of energy 
transformation systems [5] and presented in [6], the equivalence 


conditions are : i 
E{0.(t)} = E{o,(t)} 
Ca (t) = Ca (1) 
C3) =C, (T) 


where : C stands for the autocorrelation function, t is the time 
interval, and dot denotes the time derivative. The results of 
calculation are as follows : 


(22) 


O, = Op 


s .-| {scao] 


1/2 (23) 


Í o’Slodo 
o, = |7 


EO 


where : s, is the standard deviation 
of the equivalent stress amplitude. 


By analogy with the stress (6) and Eq. (13), the expected 


value of the fatigue-critical amplitude of the equivalent stress 
is as follows : 


O 


1/m 
O 
a= (0.257)? (1 - F2 —22 (24) 
ar (1 ‘i m) 
2 
where : T, is the required design life (in seconds). 
Consequently, when the following condition is met : 
(07 1/2 Oo 
Z,(1- 32) <O.5n)'s,<L(I-32) 05 
pll- Ro) < O.5n)s, R) @® 
the probability of fatigue failure during 
the time T, can be estimated as : 
2 
a 
P; =P{a 2a} exp (-=s;) (26) 
28, 
If: 
Oo 
(0.57) ?s, > L{ = =o) (27) 
Re 
the low-cycle fatigue failure may occur. 
In the case of : 
1/2 o 
0.57) s, < aN - wl (28) 
e 
an infinite fatigue life may be expected. 
Its probability is considered in [6]. 
CONCLUSIONS 


The deterministic and probabilistic fatigue-critical ampli- 
tudes of uniaxial asymmetric stresses are defined as those 
which lead to high-cycle fatigue failure during the required 
design life. 


It is found that the probability of fatigue failure under 
random-amplitude stress is equal to the probability of 
exceedance of the fatigue-critical stress amplitude. 


It is shown how the probability of fatigue failure under 
stationary random stress of known power spectral density 
can be estimated by means of the fatigue-critical amplitude 
of the equivalent random-amplitude stress. 


NOMENCLATURE 


v 2 
g 


rAS mmga 


SOU ZZS 


Hh 


o 


nw 


as 


— amplitude of the equivalent stress 

— fatigue-critical amplitude of the equivalent stress 

— autocorrelation function 

— fatigue damage accumulated during the required design life 

expected value 

— distribution function 

— imaginary unity 

— fatigue strength coefficient 

— maximum stress amplitude satisfying equation of the S-N 
curve (above which low-cycle fatigue is possible) 

— fatigue strength exponent 

— number of stress cycles to cause fatigue failure 

— required cycle number 

— probability 

— probability of exceedance of the fatigue-critical stress 
amplitude 

— probability of fatigue failure 

— tensile yield strength 

— power spectral density 

— standard deviation of the equivalent stress amplitude 


Nt 2 


TC 


agame 


a 
x wo @ @ g B® 


araaraara 


x 


e 


ex 6. 6-0. 9 


— standard deviation of the stress amplitude 
— time 
— fatigue limit under fully reversed tension-compression 


— phase angle of the stress (6) 

— gamma function 

— amplitude of the fully reversed stress 

— asymmetric stress 

— amplitude of the stress (1) 

— fatigue-critical stress amplitude 

— equivalent mean stress 

— equivalent stress 

— mean stress 

stationary random stress of zero mean value 
— amplitude of the stress (6) 

— asymmetric random-amplitude stress, asymmetric random 
stress 

time interval 

— phase angle of the equivalent stress 

— circular frequency 

— circular frequency of the equivalent stress 


* — complex conjugate 
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Application of movable approximation 
and wavelet decomposition to smoothing-out 
procedure of ship engine indicator diagrams 


Stanistaw Polanowski 
Polish Naval University 


ABSTRACT 


In this paper - on the basis of indicator diagram processing taken as an example - were 
shown possibilities of the smoothing-out and decomposing of run disturbances with the 
use of the movable multiple approximation based on the least squares criterion. The notion 
was defined of movable approximating object and constraints used to form approximation 
features. It was demonstrated that the multiple approximation can be used to decompose 
disturbances out of an analyzed run. The obtained smoothing-out results were compared 
with those obtained from full-interval approximation of runs by means of splines as well 


as wavelet decomposition with using various wavelets, Wavelet Explorer and Mathematica software. 
Smoothing-out quality was assessed by comparing runs of first derivatives which play crucial role in the 
advanced processing of indicator diagrams. 


Keywords : smoothing-out the runs, least squares method, full-interval approximation, movable approximation, 
decomposition of disturbances, cut constraints, glued constraints, riveted constraints, 
broken constraints, weighting factors, wavelet decomposition. 


INTRODUCTION 


Advanced acquisition of diagnostic information from in- 
dicator diagrams is associated mainly with determination of 
first derivative necessary for determining heat emission cha- 
racteristics. For certain purposes also determination of second 
derivative and even third one may be necessary. 

With determination of derivatives the necessity of smooth- 
ing-out indicator diagrams is associated. A. Ralston has 
defined the notion of smoothing-out in such a way that if an 
approximation maintains information on function resulting 
from measurements and fades away disturbances then it is said 
that it smooths out (levels) measurement data [8]. 

In this work - to compare effectiveness of various approxi- 
mation methods — the indicator diagram (shown in Fig. 1, 2) 
was used; it was obtained from measurements performed in 
a cylinder of Sulzer 6AL20/24 ship medium-speed engine for 
the loading parameters : n = 750 rpm, p, = 1.8 MPa. Begin- 
ning from the instant of ignition up to the end of combustion 
process, on the run of the pressure p can be observed pressure 
oscillations of large values generated by a short measuring 
channel (10 mm long) between the combustion chamber and 
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gauge membrane. The indicator diagram of p was recorded with 
0.1°owk angular resolution and 12 bit amplitude resolution. 


FULL-INTERVAL APPROXIMATION 
OF RUNS BY USING 
THE LEAST SQUARES METHOD 


Measurement data processing is aimed at assessing real 
values of measured quantities or searching for functional 
relationships whose mathematical model is either known or 
searched for. To this end, is usually applied the least squares 
method whose essence consists in the determination of mini- 

mum of the following functional : 


N 
MIN(S)=MIN ÈS- 
1 


(1) 
where : 
MIN - operator of minimum value 
N — number of elements of a measurement set 
S — sum of squares of deviations 
y; — measured values 
y;  — approximated values. 


Relations between measured, real and approximated values 
are determined by the equality : 


Yi = Yi + Eki = Vit Fa = Yi + Emit Eri (2) 
where : 
y; — real value of measured quantity 
€,, — Measuring error 
€, — approximating error 
&,, - modelling error. 


If to find an adequate mathematical model of measured 
quantity run is not possible then it is usually approximated by 
means of a linear combination of selected elementary functions. 
As an orthogonal basis power polynomials or trigonometric 
ones are most often used; the latter - in the cases where perio- 
dic runs are analyzed. Sometimes non-linear models may be 
linearized or their parameters determined by using the theory 
of experiments. 

The known difficulties in performing the full-interval 
approximation with the use of power polynomials are tried 
sometimes to be overcome by using splines. 

A drawback of such functions is their high suceptibility to 
generate oscillations appearing especially in runs of derivatives 
derived from measurement data [6] (see e.g. Fig. 1). 

As it results from comparison of the runs of the derivative 
P'eo and p’,,, the symptoms Sog and S, have been smoothed 
out for the number of nodes equal to 60. For the number of 
nodes equal to 180, the symptoms S,,, and S, appeared - but 
with significant oscillations - on the run of p’,,.. 
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Fig. 1. Comparison of smoothing-out quality of the indicator diagram p by 

means of glued 5th order polynomials and movable approximation with the 
use of the object F3 (Fig.2). Notation : Poy Pisy Psy P isg runs and their 
derivatives obtained by approximating the run p with the use of the glued 

5th order polynomial for the number of nodes given by their respective 
indices, p’, — first derivative derived from the run p by means of the mov- 
able approximation by making use of the object F3, S- — kinetic combus- 
tion symptom, S, — fuel re-injection symptom, a - crankshaft rotation angle 
[°owk] 


The oscillations are increasing along with constraint number 
increasing. As proved, no such number of nodes exists that the 
first derivative run obtained for which, could be considered 
matching the run p’, 

Worth mentioning that the smoothing-out quality of pres- 
sure run itself in the form of the runs p’,,, p’,,, can be consi- 
dered sufficient for certain purposes despite that the quality of 
determination of first and higher derivatives is insufficient. In 
the case of rather non-dynamical runs (such as e.g. those of 
pure compression, sinusoid, etc), by applying the full-interval 
approximation with splines a sufficient smoothing-out quality 
can be obtained also in the aspect of determining derivatives, 
but usually can not in the case of the indicator diagrams of 
combustion process. 


SMOOTHING-OUT THE RUNS 
AND DECOMPOSITION 
OF DISTURBANCES BY MEANS 
OF THE MOVABLE MULTIPLE 
APPROXIMATION METHOD 


An obvious way for diminishing the errors occurring at 
application of the full-interval approximation method is to 
split a given data interval into smaller ones, that leads to the 
movable approximation process. 

The movable approximation consists in determining — in 
every point n of measurement series — a value approximated 
over a movable approximation interval of a given width. The 
point n of movable approximation is called the control point. 

If the minimum value of sum of squares of deviations is as- 
sumed the approximation criterion then for P-th step (repetition) 
of approximation the criterion can be written as follows : 


i=n-k,p 
x A a2 
MIN(S,,,) = MIN > (Yeni ri) (3) 
i=n-kjp 
where : 

MIN — operator of searching for minimum 
Son — sum of squares of deviations in the point n for P-th 
. step of approximation 
Yp — values obtained from approximation for its step P 
Yo; = Y; — measured values 
Kp» Kp — parameters of the left and right end of interval 


of approximation for its step P. 


Usually the central intervals 
of approximation, i.e. kp =k, =k, are applied. 


The simplest example of the movable approximation is 
the movable average method. A natural generalization of the 
movable average method has been the algorithms of movable 
least squares approximation with the use of higher-order power 
polynomials, elaborated by Savitzky and Golay [9]. In this 
work the algorithms of movable approximation with power 
polynomials elaborated by this author [2, 3, 5], were applied; 
they were next extended by applying the approximating objects 
with constraints [4, 7]. 

In Fig.2 are shown results of the multiple approximation of 
the indicator diagram p with the use of the approximating object 
F3 being a power polynomial of 3“ order. To obtain a sufficiently 
smooth run of the first derivative p’, the approximation by using 
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160 165 170 175 180 185 190 195 200 205 210 215 220 
Fig. 2. Illustration of smoothing-out effectiveness of the indicator diagram p 
with the use of the approximating object F3. Notation : F3 — 3rd order power 

polynomial - central object, k = 16, P = 5; p’—first derivative determined 

as linear increments of the run P; Pp; pi> Pry p7 the smoothed-out run and 
first derivative after Ist step of approximation; p, , P’,;— the smoothed-out 
run and first derivative after 5th step of approximation (Pp; = Py; ps» P r3 = 
= P's ps)! Sex — kinetic combustion symptom, S, — fuel re- injection symptom. 
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the object F3 was performed five times for the same value of 
k = 16. The run of the first derivative p’ (Fig.2), determined as 
linear increments of the curve p, illustrates the scale of difficulty 
which is to be overcome to obtain a smooth derivative run not 
loaded by dynamic errors, such as the run p’ ,, (Fig.2) assumed 
in this work for the quality assessment of smoothing-out the 
run p with the use of other approximating objects and methods. 
To this end, in Fig.2 were distinguished the peaks (symptoms) 
Seg > Sp 3 On their runs conclusions as to dynamic features of 
the considered approximation method, can be based. 

The multiple approximation method can be also applied to 
decompose disturbances. The total deviation of n-th sample of 
measurement series after performance of P-th approximation, 

is equal to: 


D Ypn= Yn T Îpn (4) 


The deviation can be presented in the form of component 
deviations for distinguished steps of approximation. If the steps 
1, g, u, Pare distinguished this is equivalent to four decomposi- 
tion levels at which four runs of deviations which satisfy the 

below given equation, can be obtained : 


DYpn = Doin + DY ten + DY gun + DYupn (5) 


Other number of decomposition levels (steps) can be also 
assumed or to perform decomposition of disturbances after 
ending the approximation. Hence, auxiliary steps and main 
ones of the approximation have been distinguished. In order 
to achieve final smoothing-out to perform the auxiliary steps 
is not necessary. 

Selection of decomposition parameters as well as asses- 
sment of its results depends on expectations, i.e. on knowledge 
ofa physical character of achieved results - otherwise they may 
be of no subject-matter merit. 

The high-frequency disturbances Dp,,,, DP, of the run p 
(Fig. 2) determined in two auxiliary steps of the approximation 
with the use of the approximating objects F3a and F3b, are 
shown in Fig.3. 
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Fig. 3. The disturbances Dp p , DP p, Separated in two auxiliary steps of 
the approximation of the run p (Fig.2) with the use of the objects F3a and 
F3b. Notation : F3a — 3rd order power polynomial, k = 2; F3b — 3rd order 
power polynomial, k = 4; DP yy, = P - Prze PP pi = P -Przy Where: p p, the 
run determined from the run p with the use of the object F3a, p,,,,— the run 
determined from the run p with the use of the object F3b. 


The disturbances Dp,,,, , DP; are composed of A/D pro- 
cessing errors, disturbances within measuring lines, measuring 
gauge errors and high-frequency components involved by gas 
passages. 

The disturbances separated from the run p in five main 
steps of the approximation with the use of the object F3 are 
shown in Fig. 4. 

The runs of disturbances shown in Fig. 4 result mainly 
from interaction of the gas channel between the cylinder wor- 
king space and gauge membrane, as well as they may contain 
a part of useful signal due to inadequacy of the approximating 
object’s model. 
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Fig. 4. The disturbances Dp determined in five main steps of smoothing-out 
the indicator diagram p by using the object F3. Notation : Dp,,,, = Pp; p1- 


Prz’ DP ps4 z Prs pp -Prs pi’ DP eg H Pr ps K Pr; pp , where : Prap (Fig. 3), 
Prs po Prs ps — the runs determined in result of approximation of the run p 
y 7 by using the object F3 (Fig.2). 

The run of the deviations Dp,, can be represented by their 
three components [1]. The deviations were summed — up with 
taking into account their similarity. 

In the smoothing-out process of the run p the objects F3a and 
F3b were not used for determining the derivative p’, (Fig.2). 
However if to include the objects into the smoothing-out process 
of the run p then in the example in question the first derivative 
determined this way will be almost identical - the differences 
will not exceed +0.01% of the maximum value of p’, . Number 
of possible decompositions of deviations of measured run from 
that smoothed-out is very high for the example in question but 
such decompositions may be of no physical sense. 


MOVABLE APPROXIMATING OBJECTS 
WITH CONSTRAINTS 


The power polynomials free of constraints applied in Sa- 
vitzky-Golay filters not always ensure a sufficient quality of 
approximation. Quality of approximation and decomposition of 
runs can be formed by modifying the approximating functions 
by applying constraints onto values of left-hand and right-hand 
derivatives of the functions in the points (nodes) q, (Fig.5). 


Yau Y-(@1) Y-¢-) Yy Yı Yo Yi Yy Yei You Yu 
1 1 1 1 


K AE a W A NA A A 


I I I j I I I I 
i n 
-k Vs qx aa Ô q qv d% * k 
Fig. 5. Schematic diagram of the central approximating object. Notation : 
y, — approximating functions, {v = -u (1) u}; q,— node coordinates, 
{v = -u (1) uł; i — approximating object axis, i = 0 — control point; 
n — axis of measurement series, n =1(1)N; 
k — width parameter of approximating object interval. 


It was assumed that approximating object is central one if its 
control point is located in the mid-point of the approximating 
object. Approximating object is symmetrical if the symmetry 
of nodes, constraints and kinds of functions (coefficients) re- 
garding the mid-point of the object, takes place. 

For the sake of easiness of formulating mathematical de- 
scription of approximating object, the axis of arguments and 
that of approximating object was introduced where the control 
point coordinate is i = 0 for every successive n. Usually the 
control point is located in the mid - point of the approximation 
interval of the approximating object (Fig.5). Ifa node is placed 
in the control point then q, , but not y, , will appear in the object’s 

model, that results from the assumed indexing principle. 
Smoothing-out features of approximating object 
can be formed by setting-up its structure 
and values of its parameters including : 


=» type of approximating functions 

number of nodes and constraints 

æ kind of constraints: cut (spline), glued, riveted, broken 
ones 

æ either symmetry or asymmetry of function coefficients, 
location of control point, coordinates of nodes and their 
parameters, type of base functions. 


y 


The cut constraint takes place in a given node q, ofapproxi- 
mating object ifat least one ofthe derivatives ofapproximating 
function is cut (unconstrained) in that node : 


EAW A) Yala Eya) © 
where : 
— order of derivative 


m 
Wm— symbol of cutting (non-constraining) 
of m-th derivative. 


Such constraints are used to build splines. Obviously, 
appearance of even a single constraint of the kind in a given 
node produces an indeterminate break of function (depending 
on disturbances). 

The glued constraint in a given node q, takes place 
if the following relationships occur : 


Yin HV VGH, 
yy-1(q,) = y\(qy) 


As it results from that definition the functions to the right 
and to the left from the node q, are of different types. 


(7) 


The broken constraint consists in setting the breaking coef- 
ficient w - # | for a given m-order derivative in a given node 
q, » that is equivalent to imposing a definite discontinuity of 
the derivative in the node in question. Values of approximating 
function in nodes are assumed continuous. By using the nota- 
tion from Fig.5 the above given definition can be expressed 

as follows : 


yo GG) = Wey GG, 


Yo (Wy) = Yay) » Way #1 

The broken constraints differ from the cut ones applied 
to spline functions by that they are assumed determined ones. 

The spline constraints are free (undetermined), that results 
from cutting the selected derivatives. The simplest broken 
object is a broken line built of straight segments. 

The riveted constraints of two sections of approximating 
object’s function occurs when the overlapping of intervals of 
left and right approximating functions takes place, that can be 
expressed as follows : 


(8) 


i # la (9) 


as well as when equality of the function or one of its derivatives 
occurs even in a single node. 
In Fig.6 is shown an example of a simple riveted 
approximating object of three nodes in which values 
of the function are riveted. 
The features of the above defined approximating objects 
are illustrated by a few examples of their application 
to smoothing-out the runs [1, 7]. 


Approximating object features can be also formed by 
inserting weights (weighting functions) into approximation 
equations. The weighting methods are used in statistical data 
analysis, especially in the case of smoothing-out small data 
sets. The weights constitute a kind of constraints imposed upon 
function values. 


hl, h q2 
Fig. 6. Schematic diagram of the approximating object riveted in three 
points : y,(q,) = V441) ¥ (Uy) = V82) + Vy) =¥4,) + for the central 


43% f i 


object the following is valid: r, = -1,,17,=-1,,9;=-4)+%=9.- 


Stas [10], basing on a reference, gave an example of 
smoothing-out object which could be effective in the case 
of smoothing-out the indicator diagram. His formula for the 
smoothing-out object in question can be expressed — after 

formal transformation — as follows : 


k 
j= > wi; 

i=-k 
where the weights w, are determined from the formula : 


—1 
k {2k 2k k 
> wel 


i D2 it+k}? & 


(10) 


W. = 
i=-k 


(11) 
i+k 


where : 


k — interval width parameter 
of approximating object in question. 


As demonstrated, multiple approximation of the run p by us- 
ing the object FOw with weights has not led to any better result 
of smoothing-out as compared with that obtained by means of 

the simple movable average (object FO) (Fig.7). 
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Fig. 7. Comparison of smoothing-out quality of the run p by using the 
movable average (object F0) and that obtained by means of the movable 
weighted average (object FOw). Parameters of the smoothing-out objects: 
F0: k= 8, P = 4; FOwa: k = 8, P = 25; FOwb: k = 32, P = 7; Poar Dein 
— the smoothed-out run and its first derivative determined by means of the 
object FOWA, Pross» P rop — the smoothed-out run and its first derivative 
determined by means of the object FOwb, D — deformation . 


As seen in Fig.7, in the case of application of the object 
FOwa, as many as P= 25 smoothing-out steps are necessary 
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to obtain the result comparable with that achieved by using 
the object FO (P = 4) which is the simple movable average. 
Additionally, the deformation D appeared on the derivative 
P’ cows: The increasing of k value up to k = 32 resulted in the 
decreasing number of smoothing-out steps P = 7, but it did not 
provide any practical benefits as compared with the results of 
application of the movable average. In both the cases the deri- 
vatives significantly differ from the reference run p’, obtained 
by using the object F3. 

However the application of weighting functions changes 
smoothing-out features of approximating object and in certain 
cases their application may be justified. 


DECOMPOSITION OF RUNS BY MEANS 
OF WAVELET FILTERS 


In technical diagnostics are observed many attempts to 
apply wavelet analysis for separating — from measured runs — 
— the signals containing diagnostic information. The pressure 
run p was subjected to decomposition by using six known 
wavelet filters (offered in the Wavelet Explorer package in 
cooperation with Mathematica software) for various parameters 
of the filters [11]. 


In Fig.8 are shown the results obtained by using Daubechies 
and Shannon filters of the parameters given in their indices. 
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Fig. 8. Comparison of the results of smoothing-out the indicator diagram p 
by means of wavelet filters of two kinds with those obtained by using the ob- 
ject F3. Notation : the indices of p and p’: D — Daubechies filter, Sh - Shan- 

non filter; numerical values of the indices, e.g. 4_28, stand for parameters 

of a given filter acc. Wavelet — Explorer description; for P’ sys) 3> P 'sn200 10° 
P’ snoop 33 the index values_3; _10,; _33 stand for the values of the para- 
meter k of the smoothing-out object (3" order power polynomial) . 


As seen in Fig. 8 the first derivatives achieved from the runs 
smoothed-out by using the filters are impermissibly distorted 
as compared with the run p’,, obtained by using the object F3. 
In the case of Shannon filter the additional smoothing-out by 
means of the movable approximation method made it possible 
to obtain the run p’,, ,.. 33 Which is very similar to the run p’,,. 
In the case of the remaining filters the additional smoothing-out 
can not provide any improvement because of subsequent loss of 
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information. The first derivative deformations analogous to tho- 
se in the case of Daubechies filter, were obtained by using also 
other wavelet filters, such as e.g. : Meyer’s, Least Asymmetric, 
Spline and Coiflet filter [11]. It should be stressed that Shannon 
filter is deemed one of the worst. However in this case it was the 
only one which made it possible — after application of additional 
smoothing-out procedure — to obtain sufficiently correct results. 
The wavelet decomposition of runs provides also disturbance 
runs in the form of the so-called details. In Fig.9 are shown the 
high-frequency details Dp, „,» DPs, obtained for two first steps 
of decomposition by means of the filter Sp3_50 (Fig.8). 
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Fig. 9. Runs of the high-frequency details (deviations) DPs,» DPs,; Separa- 
ted in result of decomposition of the run p by using the filter Sp3 Si 0 (Fig.8) . 


Comparing the runs Dp, DP... (Fig.9) with the runs 
DP op Dp,,, (Fig.3) one can postulate that they are similar 
(convergent) to each other. Degree of convergence of the runs 
can be assessed by performing direct comparison, correlation 
or coherence analysis. 

In the case of low-frequency deviations (Fig.10) a notice- 
able convergence of the runs Dp, , and Dp, with the runs Dp,.,, 
(Fig.4) can be observed; however the level of their differences 
can not be considered negligible. 
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Fig.10. Runs of the low-frequency details (deviations) DP 54» PP sy 
and Dps,, separated in result of decomposition of the run p 
by using the filter Sp3_50 (Fig. 8) . 
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Analyzing the runs Dp,,,-Dp,.,, (Fig. 4) one can observed 
that in the applied display scale the waving on the runs within 
the coordinate interval up to the angle a = 170°owk, or even 
to 172°owk, are not visible. However in the case of the run 
DP., (Fig. 10) a significant waving can be observed beginning 
already from the angle a = 160°owk. The run Dp. has not any 
equivalent among the runs presented in Fig. 4. Beginning from 
the angle a = 168°owk, a significant waving is also visible on 
the run Dp, (Fig. 10), which deforms the first derivative run 
P’ 553 so t0 an impermissible degree (Fig. 8). 

Such phenomena were observed in analyzing a number of 
other runs including those much less disturbed, and in every 
case negative result was obtained of smoothing-out process for 


determining first derivative, except the case of the application 
of Shannon filter connected with final smoothing-out operation. 
The problems may serve as subjects of separate research. 


CONCLUSIONS 


O By applying the multiple movable approximation of runs 
by using objects of appropriately selected features to obtain 
a demanded degree of smoothing-out the run and to perform 
decomposition of distortions, is possible. 


O To form features of approximating objects can be used : 
cut, glued, broken or riveted constraints as well as weights. 
The features can be formed by choosing : kind of base 
functions, width of object’s interval, number and kind of 
constraints, number and location of nodes. Multiplicity 
of approximation repetitions (number of steps) influences 
smoothing-out results specially. 


O The full-interval approximation of runs by using spline 
functions may lead to significant errors in the form of ge- 
neration of waving which does not exist really. 


O The decomposition (filtration) by means of wavelets may 
be also loaded by significant errors both as to the smo- 
othed-out run and the obtained details (deviations). It does 
not contradict the usefulness of the method for e.g. image 
compression where the disturbance level involved by wa- 
velet filters may be insignificant. Applying the method for 
filtration of diagnostic signals in engineering one should be 
careful in accepting the obtained results especially if their 
confirmation by any other method is not performed. 


NOMENCLATURE 


D — run distortion, 

F3 — movable approximating object : 3" order polynomial, 
k=16,P=5 

i — axis of arguments of an approximating object 

k — width parameter of approximation interval of central object 

m  — derivative order 

n — numerical axis of measurement series, n = 1(1)N 

P — number of steps (passages) of approximation 

p — run of cylinder pressure (indicator diagram) 

p’  — first derivative of the pressure p 

Py; — first derivative of the pressure p determined with the use of 
the object F3 

q, q, — node coordinates of the approximating object, 
{v=-u(l)u} 

cx T kinetic combustion symptom 

S, ~ fe-injection symptom 

y — approximating function 

Y=y — measurement series 

y, | — approximating function, {v = -u (1) u} 


— m-th derivative of the function y, 
a — crankshaft rotation angle 
°owk — degree of crankshaft rotation angle 
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Miscellanea 


Days of Engineering 


On 23-25 November 2006, 37" Days of Engineering 
was held on the occasion of 60" Anniversary of Technical 
University of Szczecin, the oldest technical university of 
West Pomerania, as well as of many scientific technical so- 
cieties of the town, namely : that of electricians, geodesists, 
mechanical engineers and technicians, water engineers and 
technicians as well as the Federation of Scientific Technical 
Societies NOT (Naczelna Organizacja Techniczna). 


The jubilee was celebrated under the slogan : 


Youth and Engineering — a chance 
to developing the Town and Region, 


which has had a very distinct meaning. 
It was arranged by the following organizations : 


œ the Federation of Scientific Technical Societies NOT, 
Szczecin 

** Szczecin Division of the Polish Society of Electri- 
cians 

* Technical University of Szczecin 

$ The Society of Graduates from Technical University 
of Szczecin. 
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Miscellanea 


Second Edition (revised and enlarged) 
of the monograph titled : 


On Fatigue Safety of Metallic Elements 
under Static and Dynamic Loads 


written by Janusz Kolenda, Professor of Polish Naval Uni- 
versity, Gdynia, has been recently published by Naval Uni- 
versity Publishers (e-mail : b.szyniec@amw.gdynia.pl) 

The main idea behind this monograph was to formulate 
energy-based design criteria accounting for the effect of 
static loads as well as mean and residual stresses on fati- 
gue performance under either periodic or random loads. 
For this purpose the original stresses have been modelled 
by novel equivalent (reduced) stresses and Soderberg’s 

equation was used. 


Closed - form solutions consistent with Miles formula 
for the fatigue life and Parsevals formula for the power of 
periodic signals, have been obtained. The emphasis has 
been put on how to cover both the static strength conditions 
and fatigue safety requirements under combined static and 

dynamic loads. 


The book (marked ISBN 83 - 87280 - 94 - 1) 
consists of two parts divided into two 
and three chapters, respectively. 


Part I is devoted to the loading cases in which the stress 
level lies below the fatigue “safe - life” limit. 
Part II is concerned with the high - cycle fatigue regi- 
me in which a finite fatigue life is expected. 


Chapters 1 and 3 deal with design criteria 
in deterministic approach. 
In Chapters 2 and 4 probabilistic design criteria 
are developed. 
In Chapter 5 high - cycle fatigue criteria 
for anisotropic metals are considered. 


Three appendices highlight 
the following detail problems : 


. Semiprincipal stress systems 

. On the relation between probability density functions of 
maximum values of original and equivalent stresses. 

C. Application of the theory of energy transformation sys- 

tems to fatigue assessment of metallic elements under 

steady loading conditions. 


w > 


The book consists of 166 pages and includes 81 refe- 
rence sources as well as many example calculations which 
greatly enrich its didactic merits. 


ARADEMIA MARYNABKE WOJENNEJ 
See 


Miscellanea 


RINA - KORAB Award 2006 


As in the year 2005, in autumn 2006 was settled 
the competition for the most outstanding M.Sc. project 
elaborated and defended at Faculty of Ocean Engineering 
and Ship Technology, Gdansk University of Technology 

in the academic year 2005/2006. 


The competition has been arranged by : 


* the Polish Society of Naval Architects 
and Marine Engineers KORAB 
* the Royal Institution of Naval Architects - 
- RINA, UK 
* Faculty of Ocean Engineering and Ship Technology, 
Gdansk University of Technology. 


From among 9 projects submitted to the competition 
the Piotr Jatowski’s project on “Stern quarter ramp/door 
for 6700 car capacity car carrier built in Okpo Daewoo 
shipyard” won the award in question. The project was 
worked out under supervision of Prof. Zygmunt Paszota. 


Miscellanea 


75 Years of Education 


In 2006 75 years went by of education of Polish Navy 
technical corps officers, which began with establishing En- 
gineering Division in the then Polish Navy OfficerTrainees 
School in 1931. Already after a few years it bore fruits in 
the form of unaided building ships, and generally, main- 

taining high technical serviceability of naval ships. 

In present, after intensive development changes the 
education of officers of Polish Navy technical corps is 
carried out by means of modern methods at Mechanical 
- Electrical Faculty, Polish Naval University in Gdynia. 
About high quality of the education goes to show a.o. 
the fact that only in the years 1995 — 1999 as many as 
1417 technical projects and improvement proposals were 

submitted. 
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Identification of service failures of cylinder 
valves of ship piston combustion engines 


Zbigniew Korczewski 
Polish Naval University 


ABSTRACT 


INTRODUCTION 


Correct operation of self-ignition engine, which ensures 
expected performance and efficiency in steady and transient 
operational states, depends to a great extent on effectiveness 
of charge exchange in its cylinders. Quality of the process is 
demonstrated by values of coefficients of filling the cylinders 
with fresh charge, by which the so called filling efficiency is 
determined. Values of the coefficient are determined mainly 

by two factors : 


+ optimally selected distribution phases in the sense of ensu- 
ring the most favourable opening and closing angles of air 
inlet and exhaust gas ducts (cylinder filling and scavenging), 
in full range of possible engine load changes 


+ optimum velocity of air and exhaust gas flowing through 
inlet-outlet system, which ensures effective whirling the air 
flowing into cylinder. 


To ensure optimum values of the parameters during ship 
engine service it is necessary to reach the full serviceability 
state of cylinder valves which constitute its structural parts 
most thermally and mechanically loaded. Specially sensitive 
elements are outlet valves washed by exhaust gas having the 
temperature above 1000 K. In such service conditions they are 
to fulfil additional requirements as regards heat exchange and 
resistance to abrasion and impact load in high temperature, as 
well as high demands concerning corrosion resistance. 


MECHANISM OF GENERATING FAILURES 
OF CYLINDER VALVES 


During its operation the cylinder valve is forced to move 
along its spindle in the guides which undergoes friction wear. 
The process goes at high temperature of the spindle whose 


This paper presents selected diagnostic problems of charge exchange system of ship pis- 

ton combustion engines. Theoretical background of wear process of cylinder valves was 

highlighted in the aspect of identification and sources of known and identifiable states of 

unserviceability. The presented results of endoscopic examinations concern failures of 
cylinder valves of the engines installed on Polish Navy ships. 


Keywords : technical diagnostics, ship diesel engine, valve timing 


additional task is to absorb heat from valve head. As a result 
an excessive increase of radial clearance is produced between 
the guide and spindle, which leads to an undesirable skew of 
the valve resulting in loss of cylinder tightness, gas eruption, 
lubricant leakage from the spindle — guide precision pair until 
intensive wear of the entire valve unit is reached. The phenome- 
non may especially intensively develop in the case of supplying 
the engine with fuel oil of high sulphur content. The cases are 
known of completely burned-out valves as a subsequent result 
of extensive wear of valve guide [7]. In extreme case, cracking 
the valve spindle, its falling down into cylinder space and sub- 
sequent failures of the „piston — piston rings — cylinder” system 
(TPC), including piston cracking, can happen. An observable 
symptom of worn valve guides are smoked valve springs, that 
indicates a lack of tightness of combustion chamber. 

Another, often found failure of cylinder valves is a drop of 
elasticity of tightening springs and even their fatigue failures 
[4,7]. In such situation also loss of cylinder tightness can happen 
during its filling when the springs of outlet valves are cracked, 
and also during gas exhaust process when the springs of inlet 
valves are cracked. 


THEORETICAL BACKGROUND 
OF WEAR PROCESS OF VALVE GUIDES 
AND SPINDLES 


From the equilibrium condition of the forces acting on 
the valve at the initial instant of its opening (Fig.1) it results 
that the force, R, generated by distribution shaft cam pressure 
acting onto the valve spindle face, is balanced by the sum of 
the cylinder gas pressure force, P , acting onto valve head, the 
valve mechanism inertia force P, and the spring tension force 
P , according to the equation : 


R=P +P +P (1) 
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Fig. 1. Schematic diagram of the forces acting 
onto cylinder valve at the initial instant of its opening . 


Due to cam sliding over valve spindle face the additional 
friction force P, resulting in side pressure forces of the spindle 
moving inside the guide, is produced : 


P,=Rp 


where : u — sliding friction coefficient in the point A. 


(2) 


In the same time due to action of the force R in the point 
B, apart from the lateral force R,,, appears the friction force P, 
which becomes the direct cause of valve guide wearing : 


Pi = Rgh (3) 
where : u, — sliding friction coefficient in the point B. 


Values of the coefficients u and p, depend on material pro- 
perties of sliding elements, as well as on a degree of smoothness 
and dryness of contact surface. However, they do not depend 
on a size of contact surface (unless its area is so small that 

pressure force can deform it). 

As in the spindle-guide precision pair always certain radial 
clearance 6 dependent on the spindle diameter D appears, du- 
ring valve opening certain deviation, a, of coaxiality of valve 
guide and spindle is generated. However in the case of an 
excessive wear of the guide and increased radial clearance at 
the instant when cam pressures on the valve spindle, before 
the valve head becomes separated from valve seat, the gap F 
appears on one side of the seat face, that leads to loss of valve 
tightness, but on its other side the valve head is pressed onto 
the valve seat in the point C, that leads to valve spindle bending 
in the direction of action of the force P, (around the point C). 
Consequently through the developed gap hot cylinder gases 
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get out with high velocity, which results in local overheating 
the valve head and seat materials, as well as in forming local 
erosion pits in valve seat face, Fig.2. 


Fig. 2. M401A-1 engine — the outlet valve head 
of the cylinder no. 1 — traces of erosion wear of valve seat face . 


Analysing the distribution of the forces acting onto the 
valve, shown in Fig.1 one can observed that due to action of 
the friction force P, in the points B and C of the valve spindle 
- guide contact as well as the contact in the valve seat face, the 
reaction forces R, and R „respectively, appear which generate 
the bending moments M, and M, in those points. It should be 
also observed that the force system is dynamic, i.e. that in which 
load is changing continually. In quasi-stationary approach the 
force resultants can be determined by using the condition of 

equilibrium of forces and moments, as follows : 


> P,=0=> Ro,-Rg+P,=0 (4) 
hence : 
Re =R, -P, (5) 


$ P, =0>Rc,+P, +P, +P, -P -R=0 (6) 


Assuming that at a given instant the sum of the forces acting 
onto the system is constant one obtains the following : 


P =P tP tP R 
Then the sum of bending moments respective 
to the point C is as follows : 


S Mec=0>P,(2+8)+ 


2 
-P,(S'+2+8)+RgL—-PL, =0 


(7) 


(8) 


hence : 


D D,+D+26 
p, (28) Ry 9 )-L|-RuL,=0 
2 (9) 
The reaction force in the point B, which decides on the 
rate of wear of the upper part of the valve guide, is then the 
following : 


_P (D +28)-2RLp 
B u (D+D+28+2L) 


Making use of the relation (7) one finally 
obtains the following : 


p — Pet Pot PAD, +28) —R(D, +25 + 2L,p) 
B u,(D, +D+28+2L) 


R (10) 


(11) 


When the relation (3) is taken into consideration one can 
determine the value of friction force in the upper part of the 
valve guide : 


p — Pet Pot P, (D+ 25) -R(D, + 25 + 2L) 
a D,+D+26+2L (12) 


Making use of the relation (5) one can also determine the 
reaction force in the point C, which decides on the rate of wear 
of the valve seat face : 


Rœ =R,- Rp (13) 


Considering the system ofthe forces acting onto the valve at 
the instant of separation of the valve head from the valve seat, 
Fig.3, one can observe that the spindle presses the lower part of 
the valve guide in the point D (the support point C disappears). 


Zones of the greatest 
wear of valve guide 


Fig. 3. Schematic diagram of the forces acting onto the cylinder valve 
at the instant of its full separation from the valve seat . 


Under influence of the force R in the points B and D, apart from 
the transverse forces R, and R, the respective friction forces 
P, P, are induced, which become the direct cause of the valve 
guide wearing in its upper and lower part : 


Pa = Rg, (14) 
Po = Rph, (15) 
where : 


u, — sliding friction coefficient in the points B and D. 


In such situation the reaction forces R, and R, resulting 
from the friction force P, generate the respective bending 
moments M, and M, in the points B and D. Making use of the 
conditions of equilibrium of forces and moments one obtains 
the following : 


> P,=0>Rp-Rg+P,=0 (16) 
hence : 
R,=R,- Ru (17) 


> P= 0=>-P, +P, +P, +P,-P,.-R=0(18) 


By assuming that at the given instant the sum of the forces 
acting onto the system is constant the following is obtained : 
Eo as tP tP -R (19) 
On insertion of (19) to the relation (18) 
the following is yielded : 
Pa oly ls 0 (20) 


When the relations (14), (15) and (17) are taken into account 
the reaction force in the point B, which decides on the rate of 
wear of the upper part of the valve guide : 


Pw + Ruy H 
2u 


as well as the reaction force in the point D, which determines 
the rate of wear of the lower part of the valve guide, can be 
determined : 


- Pa — Ruu 
Ro= 2u 


In Fig. 3 the zones of the greatest rate of wear 
of the valve guide are marked red. 


Rg= (21) 


(22) 


Summing up the above performed considerations one can 
enumerate the design factors determining the reaction forces R, 
R, and R,,, which detrimentally influence the design structure 
of valve spindle guides : 


> the height of spindle above valve guide, 1, which should be 
as small as possible 

> the valve guide length L, as well as the distance from the 
upper edge of the guide to the seat face of valve head, L, 
which should be as big as possible 

> value of the friction force P generating side pressure forces 
of the spindle moving along the guide, which should be kept 
as small as possible. 


The last design factor decreasing wear of cylinder valves 
is controlled by changing the character of valve rocker friction 
against the valve spindle face — from sliding friction to rolling 
one. As shown in Fig.4 it is possible by mounting the rollers at 
the rocker’s end cooperating with the valve spindle. 


Fig. 4. Schematic diagram of the forces acting in the “ valve rocker 
with roller — face of valve spindle” system . 


Introducing simplifications, one can assume that the roller 
which substitutes the cam, generates only the pressure force 
R applied to the valve, in absence of any friction force P (its 
value is about ten times smaller than that in the case of sliding 
friction) [5]. Then no side components responsible for wearing 
the guide will be produced. The system will be in the state of 
equilibrium described by the relation (1). However because of 
the drop of elasticity (relaxation) of the valve spring, its service 
wear, wear of spindle face and valve head, certain transverse 
reaction leading to guide wearing will be always produced and 
the system will reach the equilibrium shown in Fig.1, except 
that in the point A no friction force P, will be present. 
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From the condition of equilibrium of forces and moments 
the reaction forces are determined as follows : 


Reo =R, (23) 
Ra tP, tP, +P -P - R50 (24) 
The friction force in the point B, which determines valve 
guide wearing, is described as follows : 
(25) 
where : u, — sliding friction coefficient in the point B. 


Assuming as before that at a given instant the sum of forces 
acting onto the system is constant (7) one obtains the expres- 
sion which describes the sum of bending moments respective 
to the point C : 


(26) 
_ ea R,L=0 
hence : 
Py(oi+8)-Ry m (2P) 1} 0 (27) 


Hence the reaction force in the point B, which decides on the 
rate of wear of the upper part of valve guide, is as follows : 


_ P D42) 
B wD +D+28)-2L 


It is quantitatively equal to the reaction force in the point C, 
which decides upon the rate of wear of valve seat face. 


=Ro, (28) 


The expression which describes the sum of bending mo- 
ments respective to the point B, obtains the following form : 


>My- oR T 


(29) 


D+26 
=0 
a 
Hence the vertical component of the reaction force 
in the point C is expressed by the following formula : 
R — Rex —Py(D +28) 
Cy D,+D+26 

Applying the analogical considerations to the system shown 

in Fig. 3 where the roller is used instead the cam to neglect the 


friction force P, in the point A, one can determine the equilib- 
rium equations for the following forces : 


> Fs=0SRa= Ry 


+R,L-P,,( 


(30) 


(31) 
> Py=O=>-P, +P ,+P,+P,-P,-R=0 (2) 


Assuming that at a given instant the sum of forces acting 
onto the system is constant one obtains the following : 


Pa Pi ie ak (33) 
Inserting it into (32) one obtains : 
Po = b= 2 = 0 (34) 


After taking into account the conditions : 


Pi = WR, Es = WR, = Pi = Po (35) 


and 
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one can determine the reaction forces in the points B and D, 
which decide upon the rate of wear of the upper and lower part 
of valve guide : 


Pa _ Pe Pa _ Pa + Ph +P, -R 
Ho Mm 2i 2u; 


From comparison of the expressions (21), (22) and (36) 
which describe the reaction forces in the support points B 
and D, for the system fitted with the distribution shaft cam 
and that with the valve arm roller, respectively, it results that 
the use of the roller makes the rates of wear of upper and 
lower part of valve guide equal. And, the greater the radius 
of the roller the smaller the rate of wear. At the initial instant 
of valve opening the rate of wear of valve seat face will be 
also lower due to a significant decrease of the reaction force 
in the point C. 


EXAMINATIONS OF CYLINDER VALVES 
OF SHIP DIESEL ENGINES 
DURING SERVICE 


R= 


=Rp 66) 


Novel methods of diagnostic tests are commonly introduced 
into operation process of ship diesel engines. The dynamically 
developing endoscopy which has been earlier applied only 
in medical examinations, now serves as a very useful, even 
indispensable tool for assessing the technical state of complex 
ship engines. 

Endoscopy is a disassembling-free method for visual - opti- 
cal inspection of interior of machines and devices by means of 
specular instruments (endoscopes). 

For endoscopic examinations of the engines installed 
on Polish Navy ships the following instruments are used : 
IF8D4—15 fiberscope and a kit of borescopes of OLYMPUS 
and STORZ firms, differing to each other by the length of 
optical system, its diameter and angle of observation of a diag- 
nosed element, namely : 90cm/8mm/90°, 55cm/8mm/90°, 
45cem/8mm/90°, 50cm/6mm/90°, 30em/4mm/0°, 30cm/10mm/ 
120°- see Fig.5. The instruments make it possible to examine 
and prepare photographic documentation of engine’s internal 
elements, through inspection openings having their diameter 
greater than 5 mm. A special digital photo-camera, Camedia 
C-—2500L made by OLYMPUS firm, is used to perform dimen- 
sional analysis of detected failures, visualize them and record 
in a data base. 


Fig. 5. The endoscopic diagnostic system of OLYMPUS firm : 
1 — kit of borescopes, 2 — fiberscope, 3 — kit of light sources, 
4 — digital photo-camera, 5 — photo-printer . 


The camera is connected with a borescope or fiberscope 
by means of special links (adapters). 


The length of the fiberscope’s elastic light pipe whose 
controllable end makes observation in an arbitrary direction 
possible, is equal to 1500 mm. It has replaceable ends making 
observation within front and side sectors of different obser- 
vation angles possible. Owing to this, to a great extent are 
increased manual possibilities of inspection of interior of air 
and exhaust gas flow passages of engine and turbo-compressor 
system. 

Borescopes of different lengths and rigid optical system 
make it possible to carry out observations within front and 
side sectors in a wide range of variability of observation angle. 
The 30cm/10mm/120°optical system is especially useful in 
diagnosing the engine combustion chambers, the valve seats 
fixed in lower plate of engine head in particular. Borescopes are 
also very useful during inspection of guide vanes and moving 
blades of turbo-blower. 

In Fig.6 is presented a way of conducting the endosco- 
pic examinations of ship engine cylinder systems by using 
a borescope and fiberscope. And, in Fig. 6 is presented a way 
of getting access to interior of the cylinder liner of the ship 
diesel engines : M401A-1(2) and 16V149TI Detroit Diesel, 
for endoscopic examinations. 


a) 


Fig. 6. Ways of inserting the end of borescope and fiberscope into cylinder 
interior of : a) M401A-1 diesel engine — through the holes remaining from 
dismounted injectors, b) 16V149TI Detroit Diesel engine — through inlet air 
windows in cylinder liner . 


When the injector is dismounted the borescope (fibersco- 
pe) makes it possible to assess technical state of piston head, 
cylinder liner surface, cylinder head and other units fixed in it 
such as spray nozzles of remaining injectors, inlet and outlet 
valves, starting valves etc (Fig.7). The endoscopic method of 
inspection is especially useful in diagnosing multi-block and 
multi-cylinder engines. For instance, in the case of radial en- 
gines, e.g. the M503 or M520 engine already installed in the 
engine room, an access to its lower mono-blocks and lower 
parts of reduction — reversing gear is very difficult. In order 
to perform inspection of the engines together with their gears 
they must be uncoupled from propeller shaft, next inclined and 
lifted and sometimes even rotated inside the engine room so 
as to make the first or seventh cylinder block accessible. From 
operational practice it results that a fiberscope of a sufficiently 
long light pipe makes it possible to avoid those inconveniences 
and thus to save execution time of overhauls and associated 
costs even by 25-30% [7]. 


a) Access to cylinder liner interior b) 
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Fig. 7. Endoscopic examination of ship engines — access to cylinder liner 
interior: a) 16V149TI Detroit Diesel engine, b) M401A-1(2) diesel engine . 


FAILURES OF CYLINDER VALVES 
OF SHIP ENGINES 


Systematic endoscopic examinations of ship engines 
installed on Polish Navy ships 
are carried out in the following situations : 


** during preventive surveys (once a year at least) 

œ at current assessment of technical state of an engine when 
to prolong the time between repairs is necessary 

% in the case of an increased vibration level, found metal 
particles in lubricating oil, occurrence of sudden deviations 
from the trend line of mean indicated pressure (indicated 
power) in a cylinder, increased exhaust gas temperature, 
excessive smoke emission etc. 

% when to dismount the engine head is difficult and time- 

-consuming. 


On the basis of multi-year endoscopic examinations of 
ship engines a procedure of their technical state assessment 
in service conditions was elaborated. It contains the necessary 
scope and scheduling principles of examinations of engine 
interior, aimed at detection of possible defects of particular 
elements of functional systems of the engine. For every type 
of the engines used on Polish Navy ships detail instructions for 
realization of diagnostic examinations by means of the fiber- 
scope and kit of borescopes, were elaborated. The identified 
failures are photographically recorded to make documentation 
of the detected defects and to determine a tendency of their 
development. Results of the examinations are saved in a com- 
puter data base [1,2,3,7]. 

High effectiveness of the method in question, at a rela- 
tively simple use of the diagnosing instruments, has been 
confirmed during almost 15 years of experience in endoscopic 
examining the ship engines, gained by specialists from Polish 
Naval University. As a result of the performed examinations 
were identified many material defects which could be trans- 
formed into serious hazards to engine reliability in the case of 
their further development. 

The selected failures of cylinder valves of ship diesel en- 
gines, identified during endoscopic inspections in service, are 
presented in Fig.8 and 9 [1,2,3 and 7]. 

The heavy oiling-up of inlet and outlet valves goes to show 
that the valve spindle guides are excessively worn. In order not 
to allow for premature wearing-out the guides and for bending 
the valve spindles, that always leads to loss of tightness of 
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a) Inlet valve head in the cylinder head of 5BAH22 
engine — deposit layer on valve head surface 


b) Outlet valve of 16V149TI Detroit Diesel engine — 
— lubricating oil and tar compounds deposited on valve head 


c) Valve head of opened inlet valve in the cylinder head seat 
of 75H12 WOLA engine — smooth valve seat face, without 
any traces of wear (reference state) 


d) Valve head of opened outlet valve in the cylinder head seat 
of 75H12 WOLA engine taken out of service — characteristic 
wear sill on valve seat face 


e) Outlet valve of M503A engine — loss of valve tightness 
due to bent valve spindle 


J) Inlet valve head in the face of cylinder head of M401A-2 
engine — oil layer (visible oil drop) on valve head surface, smooth 
valve seat face, without any traces of wear and corrosion 


Fig. 8. Defects of ship diesel engines, identified during endoscopic examinations . 


cylinders, two values of radial clearance between the guide 
and spindle are provided : 


X optimum one, i.e. that for assembling 


ye ultimate one, i.e. that qualifying a given valve for replace- 
ment or regeneration. 
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It is recommended to keep the assembling values of radial 
clearance of inlet valves, related to the valve spindle diameter 
D, within the range of (0.004 + 0.006) D and their ultimate 
values — within the range of (0.01 + 0.015) D, which means 
that the radial clearance value may increase even 2 + 3 times 
during engine service. In the case of outlet valves the respective 


values of radial clearance amount to : (0.006 + 0.01)D — for valve guide. A measure of the bend can be the angle a which 
assembling, and (0.015 + 0.025) D — ultimate ones [6]. determines coaxiality of the spindle and guide, or the gap F 

The ultimate value of radial clearance results from the which appears on the valve seat face. (Fig. 1 and 3). As results 
permissible bend (curvature) of valve spindle moving within from the system of forces acting onto cylinder valve, shown in 


a) Outlet valve head in the cylinder head of 5BAH22 engine — b) Inlet valve head in the cylinder head of 5BAH22 
— deposit layer on valve head surface engine — smooth valve head surfaces, without any traces 
of contamination, visible surface roughness on valve seat face, 
resulting in exhaust gas blow-by 


c) Outlet valve head in the cylinder head of Wartsila engine — d) Outlet valve seat of Wartsila engine — carbon deposit layer close 
— carbon deposit layer close to valve seat face to valve seat face — exhaust gas blow-by through non-tight valve 


e) Outlet valve head of M401A4-1 engine — oily deposit layer J) Outlet valve head of M401A-1 engine- carbon deposit 
on valve head surface and tar compounds on valve seat duct surface, valve seat 
face with some traces of carbon deposit 


Fig. 9. Defects of ship diesel engines, identified during endoscopic examinations . 
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Fig. 3, the angle of valve spindle bend is equal to the angle of 
skew of valve head in compliance with the expression: 


a-25 F 


Lp D; 


(37) 


The quantities are determined, at the same values of radial 


clearance, by the following design parameters : 


=> guide length 

= valve spindle length 

= valve head diameter 

= mutual location of the spindle and guide 


From the expression (37) it results that the dimensions of 
the gap of valve seat face, F, increase when the radial clearance 
between guide and spindle, 6, increases, the greater the valve 


head diameter D,, and the smaller the guide length L, 


25D, 
Lp 
From the relation (38) it also results that in the case of the 
predetermined permissible value of the valve seat face gap, 
F, the greater the guide length L and the smaller the valve 


head diameter D, the greater the ‘ultimate value of the radial 
clearance ô : 
_ FLp 


B= 2D, 
Additionally, worth mentioning that the ultimate value 


of the radial clearance does not depend 
on the valve spindle diameter D. 


F= (38) 


õ 


(39) 


For the considered ship diesel engines the radial clearance 
values of inlet and outlet valves are presented in Table [6]. 


Table. Assembling and ultimate values of the radial clearance between the spindle 
and guide of inlet and outlet valves of the engines : M401A-1(2) and 16V149TI Detroit Diesel . 


Inlet valves 


Outlet valves 


Type Assembling 


of Valve spindle 
engine diameter 
[mm] 


clearance 
value 
[mm] 


M401A-1(2) 18.00 0.094 


Ultimate Valve Assembling | Ultimate 
clearance spindle clearance clearance 
value diameter value 
[mm] [mm] [mm] 


0.144 


16V149TI Detroit 
Diesel 


Gap-valve distribution system 


CONCLUSIONS 


O Leakage of lubricating oil from valve guide is a symptom 
of dropping compression pressure in engine cylinder as 


a result of loss of tightness of combustion chamber. 


O Rate of wear of cylinder valve guide depends on pressure 


force exerted by valve spindle onto valve guide wall, 


on their mutual sliding. 


and 


O And, velocity of the sliding is directly proportional to the 


mean piston velocity and crankshaft rotational speed. 


NOMENCLATURE 

D — diameter of valve spindle 

D, — diameter of valve head 

F — height of valve face gap 

LE — height, length 

M — moment of force 

P, R — force 

p — pressure 

a — angle 

ô — radial clearance between valve spindle and guide 
u — sliding friction coefficient in the point A 

H — sliding friction coefficient in the points B, and D 
T — time 


Indices and abbreviations 


— of inertia 

— of gases 

— an ultimate quantity 
— of guide 

— of spring tension 

— of friction 

— aresultant quantity 


er espns 
4 
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X,Y — a component quantity 

Z — of valve 

A, B, C, D — points of contact (support) of valve system elements 
TPC — the,, piston — piston rings- cylinder” system 

ZS — self ignition 
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ABSTRACT 


Service speed obtainable by a ship in real weather conditions when sailing on a given shipping route, is 
one of the major parameters which have great impact on ship operation costs. The so far used, very ap- 
proximate method of service speed prediction based on “service margin”, is very little exact. In this paper 
a new method based on additional ship resistance dependent on mean statistical parameters of wave and 
wind occurring on a given shipping route, is presented. The mean long-term service speed is calculated on 
the basis of the calculated additional resistance and the screw propeller and propulsion engine parameters. 
Also, a new definition of service margin and a way of its calculation is presented apart from the results of 
the mean service speed calculation depending on ship s type and size and shipping route. 


Keywords : ship service speed, wind, waving, shipping route, service margin, long-term prediction. 


PARAMETERS OF THE SHIPS 
FOR WHICH CALCULATIONS 
WERE PERFORMED 


Calculations of the mean additional resistance and mean 
long-term service speed were performed for 7 existing ships 
of determined propulsion systems and 15% service margin. 
They represented two types of ships : bulk carriers and conta- 
inerships of different size. The main parameters of the ships 
are given in Tab.4. 

All calculations for each ship were performed for con- 
stant ship’s draught and constant state of hull and propeller 
surface, that made it possible to more accurately compare the 
ships sailing on a given shipping route, by taking into account 
ship’s type and size. The used computer software has been 


so designed as to make it possible to perform calculations at 
variable draught (considered as a random variable of a given 
occurrence probability) and for changeable state of hull and 
propeller surface. 


SHIPPING ROUTES 


For the calculations of the mean additional resistance and 
mean long-term service speed the most representative shipping 
routes were selected on the basis of [1]. They are specified in 
Tab.5 and shown in Fig.16. 

All additional data which describe ship’s voyage on a given 
shipping route (i. e. specification of sea areas crossed by the 
route, probability of ship’s staying in a given sea area, f,, ship’s 
course angle y and probabilities of ship’s courses in a given 
sea area, f) can be found in [4]. 


Tab. 4. Ship main parameters . 


Container- | Container- | Container- Bulk Bulk Bulk Bulk 
-ship -ship -ship carrier | carrier | carrier | carrier 
K1 K3 M1 M2 M3 M4 
Length b.p. 
Breadth B [m] 22.30 25.30 32.24 23.0 25.3 32.2 32.2 
Draught T [m] 8.25 9.85 10.50 8.55 10.65 12.05 11.6 
Displacement at T V [m°] 17300 29900 47250 21441 40831 56396 73910 
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> 
S| | Contract speed | V, [m/s] 8.44 9.62 10.80 715 | 7.72 8.2 8.28 
=| | Propeller diameter | D, [m] 52 6.15 742 5.0 5.8 6.2 6.5 
= Propeller pitch P [m] 3.9 5.8 6.8 4.2 3.9 5:7 4.5 
S| Smmm N kw] 6930 13750 26270 4710 | 7500 | 12000 | 12720 
power n 
Digawe TA 2.33 1.82 1.72 185 | 203 | 164 | 193 
rated speed 2 
Still - water 
resistance of ship | R [kN] 461.3 858.8 14775 | 4059 | 600.4 | 865.6 | 947.6 
for T and V, 
Service margin 
assumed in 
e jeg 15 15 15 15 15 15 15 
design 


CALCULATION RESULTS 


The calculation results for every ship and shipping route are 
presented in the form of bar charts for the additional resistance 
AR and service speed V,,. The speed the ship is able to achieve, 
was calculated under the assumption that its set value has to be 
maintained or its maximum value achievable in a given weather 
conditions at simultaneous keeping the set ship course, has to be 
determined. The obtained bar charts were also approximated by 
using the distribution functions : that of the additional resistance 
f(AR), and that of service speed f(V,,), with the determined R? 

(share of the “resolved” variance), [2]. 


Under the bar charts for a given ship 
and shipping route the following data are tabularized : 


> the maximum value of additional resistance, AR ax Which 


occurred for a given ship on a given shipping route 
> the mean value of additional resistance, AR, calculated from 
the formula : 
NR 
>; Prr; (AR, = const) 
AR= += (45) 
Nr 


> the relative resistance increase, PR, calculated from the 


In each case the calculations were performed for both legs formula : 
of voyage (there and back) hence every shipping route has two AR 
indices : a and b (e.g. shipping route no. 2 — 2a and 2b). PR= R 100% (46) 


Tab. 5. Specification of the shipping routes used for calculation . i 
of additional resistance and service speed of ship . > the set value of ship service speed, Vig 


Nannie > the minimum value of ship speed, V» which occurred on 


: a given shipping route at keeping a set course 
1 South America — West Europe > the mean long-term value of ship speed, Vg , on a given 
2 US East coast — West Europe shipping route, calculated from the formula : 


US East coast — Gulf of Mexico — West Europe ny 
US East coast — Mediterranean — West Europe E > Pry, Vi (AR; = const) 
Indonesia — Japan VE= H = (47) 
Persian Gulf — Japan > Pry, 
North Africa — West Europe i= 
North Africa — US East coast > the probability, Py, , of maintaining the service speed, V, , 
Persian Gulf — Africa — West Europe on a given shipping route at keeping a set course 


> percentage of the events in which ship speed has been redu- 
ced for the reason of exceeded ship’s sea-keeping criteria 
percentage of the events in which no result of calculations, 
West Europe — South and Central America i.e. no value of ship speed, has been found. The lack of solu- 


West Europe — Mediterranean — Persian Gulf — Japan 
West Europe — Panama Canal — US West coast > 
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tion resulted from the fact that it was assumed to keep a set 
course, that was not possible in very high and oblique waves 
- the ship course should be changed in such situations. 


The calculated bar charts for all the ships (Tab.4) and all 
the shipping routes (Tab.5) together with the remaining results 
can be found in [4]. The example bar charts for K1 ship (Tab.4) 
sailing there and back on the chosen route (South America — 

— West Europe), are presented in Fig.17. 


In Fig.19 are showed the values of additional resistance 
index for the same ship, defined by means of the formula 
(46) depending on a shipping route. The index should not be 
identified,as for its value, with the service margin because 
it constitutes a definite long-term statistical quantity tightly 
connected with parameters of weather conditions occurring on 
a given shipping route, hence value of the index is changeable 
whereas value of the service margin assumed in design process 
is of a constant value. However, value of the index may suggest 


Ship : K1 
Shipping route no. la 


Histogram and distribution function 
of additional resistance 


y = exp[(-0.5831113) + (-0.024341)x] 
R’= 0.986 


Histogram and distribution function 
of service speed 


y = exp[(-64.2556) + (7.60567)x] 
R’= 0.999 


0.6 
0.4 
0.2 
o0 mM N 0 
ron EA EA 2S © a + a <+ S Re} a o0 + 
Ss 3 & f $ g E g ao 6 6 E n & 
— al N N fog) 

Still-water resistance at the speed V, [KN] R | 461.26 Set service speed [m/s] Vi, | 8.44 
Max. additional resistance [kN] AR ax | 367.70 Minimum speed [m/s] V ain | 1-88 
Mean additional resistance [kN] AR | 35.77 Mean speed [m/s] Va 8.38 

Resistance increase [%] 7.8 Probability of V, Py 0.90 
Percentage of no-solution events [%] 0.11 Percentage of speed-reduction events [%] 0.11 
Shipping route no. 1b 

0.7 1 
06 y= exp[(0.1673814) + (-0.055117)x] y = exp|[(-73.50207) + (8.703594)x] 
R’= 0.994 0.8 R’= 0.999 
0.6 
0.4 
0.2 
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Still-water resistance at the speed V, [kN] R | 461.26 Set service speed [m/s] Vi, | 8.44 
Max. additional resistance [kN] AR nax | 365.74 Minimum speed [m/s] V nin | 1-88 
Mean additional resistance [kN] AR | 28.02 Mean speed [m/s] V; 8.37 

Resistance increase [%] 6.1 Probability of V,,. Pe 0.92 
Percentage of no-solution events [%] 0.22 Percentage of speed-reduction events [%] 0.07 


Fig. 17. Bar charts and distribution functions of additional resistance and service speed for K1 ship . 


ANALYSIS OF THE ACHIEVED RESULTS 


The presented calculation results clearly show which value 
of the mean long-term service speed can be achieved by the ship 
or which value of probability of maintaining a given service 
speed on a given shipping route can be achieved at a service 
margin value assumed in design process of ship propulsion 
system (e.g. 15 %). Both the quantities for K1 ship on different 
shipping routes are shown in Fig.18. 


how large service margin for a given shipping route should 
be to maintain the set service speed at a given probability of 
its exceeding; it can be clearly observed when comparing the 
diagrams in Fig.18a and Fig.19. It means that in the case of 
applying the traditional method of determining the propulsion 
engine output power, with the use of service margin, percent- 
age value of service margin for a given ship should be tightly 
correlated with statistical parameters of weather conditions 
on a given shipping route, and with an assumed probability of 
maintaining an assumed long-term service speed of ship. 
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Fig. 18. Ship : K1. 
a) probability of maintaining the set service speed V, 
b) mean long-term service speed . 


SS aa 


Fig. 19. Index of additional resistance increase for K1 ship, 
depending on shipping route . 


In Fig. 18 and 19 is presented an influence of shipping routes 
on the parameters associated with speed of K1 ship. The same 
relations for the remaining ships can be found in [4]. 

In the next figures an influence of type and size of ship on 
the parameters associated with its speed on a given shipping 
route, is presented. As the considered ships differed to each 
other with regard to the contract speed Ve the following index 
of maintaining the set speed was introduced : 


wea V 
ZE 

The index makes it possible to more accurately compare 
performance of particular ships on a given shipping route. In 
Fig. 20 + 25 are presented : the probability, Py 3 of maintain- 
ing the set service speed V,,,., the index WP acc. (48), as well 
as the index of increase of additional resistance acc. (46), for 

bulk carriers and containerships, respectively. 


100% (48) 


Fig. 20. Probability of maintaining the set service speed V ,, 
of bulk carriers : blue - M1; red - M2; yellow - M3; green - M4. 


30 POLISH MARITIME RESEARCH, No 2/2007 


Fig. 21. Index of maintaining the set service speed V, of bulk carriers : 
blue - M1; red - M2; yellow - M3; green - M4. 
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Fig. 22. Index of additional resistance increase for bulk carriers : 
blue - M1; red - M2; yellow - M3; green - M4 . 


Fig. 23. Probability of maintaining the set service speed V, 
of containerships : blue - K1; red - K2; yellow - K3 . 
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Fig. 24. Index of maintaining the set service speed V, 
of containerships : blue - K1; red - K2; yellow - K3 . 
From comparison of the diagrams it results that the perfor- 
mance of the ships ona considered shipping route differs to each 
other. Some ships show a greater value of the probability, Pyp 


of maintaining the set service speed V, and a greater value of 


ZE? 


the index WP than those showed by other ships. It speaks for 
different effectiveness of their propulsion systems. 


Fig. 25. Index of additional resistance increase for containerships : 
blue - K1; red - K2; yellow - K3 . 


In the traditional designing of ship propulsion system and 
determining of service margin value (e.g. on the level of 15% 
[3,4]) the set service speed in real weather conditions is as- 
sumed lower than the contract speed by about 1 knot (for ships 
of the contract speed up to 20 knots). Hence for K1 ship were 
performed calculations of the probability, P, of maintaining 

a set service speed for its two values : 


aVV 
2V,,=V,— 1 [kn] 


and their results presented in Fig. 26. 


Fig. 26. Probability of maintaining the set service speed V,, of K1 ship . 


In the second case (i.e. V, = V, —1 [ kn]) the probability 
Py, for many shipping routes is close to 1, that is in compliance 
with the current designing practice and to a certain extent con- 
firms (and verifies) accuracy of the performed calculations. 


A PROPOSAL OF DEFINITION 
OF SERVICE MARGIN 


On the basis of the performed calculations and conclusions 
resulting from the diagrams in Fig. 18 + 25 was elaborated 
a proposal of definition of the service margin which associates 
a set service speed with probability of its maintaining on a given 

shipping route during a long time period. 


Proposed definition of the service margin k, 
(which determines surplus 
of engine output power) 


For a given ship sailing on a given shipping route, the ser- 
vice margin k, should have such percentage value as the ship 


in assumed loading conditions and assumed state of its hull 

and propeller surface, would be able to maintain the service 

speed, V,,., assumed in the frame of long-term prediction, at 
the assumed exceeding probability P 


The definition concerns a designed ship intended for sailing 
on a given shipping route, however it can be also generalized 
to cover an arbitrary sailing region. 


It should be observed that : 


% too large value of V,, relative to V, (contract speed) or too 
large value of the probability P „ may cause that k, would 
take large values and such ship would appear operationally 
unprofitable for its owner 

% areasonable (smaller) value of V,,, may cause that at a rea- 

sonable value of k, , the probability P,,, would reach the 

value of one (P; = 1). 


The above given comments mean that ship owner, 
when ordering a ship, should be aware of : 


# a value of the long-term service speed he expects to be 
reached and 

Æ at which level of its maintaining probability the designed 
ship would appear operationally profitable for him. 


For the so defined service margin were performed calcula- 
tions aimed at determination of its percentage value depending 
on shipping route under assumption that the set service speed 
V,,» maintained during long time period, will be equal to V, 
(Vz = Vo with the probability P,,. = 0.95. 


The results of service margin calculations for K1 
ship in function of shipping route are presented in Fig.27. 
407 k; [%] 


Fig. 27. The service margin k, calculated for K1 ship in function of shipping 
route, for V, = V and P, = 0.95; blue - service margin for one-leg voyage 
of the ship; red - mean service margin for a shipping route . 


CONCLUSIONS 


From the performed calculations 
the following conclusions can be drawn : 


O Additional resistance of ship and its mean service speed 
depend on statistical weather parameters of shipping routes, 
hence ship service parameters depend on a kind of shipping 
route hence full utilization of statistical long-term parame- 
ters of waves and winds makes more accurate predicting 
ship service speed on a given route, possible. 


O Influence of ship’s type and size on its service parameters 
(resistance, speed) is observable but not so great as that of 
shipping route, but definite differences in effectiveness of 
the same ship on a given shipping route can be observed, 
that speaks for different quality (perfection) of designed 
propulsion systems. 
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Service margin’s value can be determined on the basis of 
ship service speed assumed by ship owner and probability 
of its maintaining on a given shipping route, hence it is 
not to have a constant value for a given shipping route, but 
it may result from targets and strategy which ship owner 
would apply to a given ship. 


On the basis of the presented method was elaborated 
a computer software which may be also used 
for the following purposes : 


Comparison and assessment of ship designs as well as exi- 
sting ships from the point of view of their service costs 


Optimization of a shipping route for a ship, especially in 
rough weather conditions 


Investigation of many other service parameters of ship and 
its equipment, e.g. frequency of using ship’s steering gear 
on a given shipping route. 


Calculation results obtained in the frame of long-term 


prediction of ship service parameters will be presented in 


separate papers. 


NOMENCLATURE 


f(V,) — probability distribution function of occurrence of the 


instantaneous ship speed V; 


f(AR,)— probability distribution function of occurrence of the 


f 


A 


f 
v 


instantaneous additional resistance AR, 

— probability of ship’s staying in a given sea area A 

— probability of ship’s sailing with the speed V and course 
angle y 

— service margin 

— number of intervals for additional resistance 

— number of intervals for ship speed 

— relative increase of resistance 

— total occurrence probability of the additional resistance 
AR of a given value 


P~ — total probability of developing the speed V by a ship at 
occurrence of the additional resistance AR 


P, ~ probability of developing by a ship the set service speed V 

R — still-water ship resistance 

R? — share of the “resolved” variance 

V z Service speed 

V, — average long-term service speed of a ship sailing on 
a given shipping route 

V, — contract speed 

Vg 7 Set service speed 

WP -— index of maintaining the set service speed 

AR — additional resistance resulting from rough weather 
conditions 

AR ~ statistical average value of additional resistance 

y — geographical angle of ship course. 
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The Foundation for Safety 
and Environment Protection 


25 years of training on manned ship models 


In 1981 in Itawa (Silm lake) on the basis of agree- 
ment between the then Maritime High School, Gdynia, 
and Gdansk University of Technology was commenced 
regular education activity aiming at training ship masters, 
first officers and pilots in ship manoeuvring. The trainings 
were based on large manned floating ship models equipped 
with devices simulating operation of the ship systems 
engaged in manoeuvring. The first training model rep- 
resenting a tanker is one out of nine models used today 
for training courses. As the years go construction and 
outfitting of models as well as training water areas, at first 
simple and sparse, have been becoming more and more 
perfect and better adjusted to fulfilling their role. Also, 
organizational form of the training has been changed - in 
1989 the Foundation for Safety and Environment Protec- 
tion, Iława, was established by both the universities with 


the aim of carrying out trainings and research in the area 
of safety at sea. 


In spite of many difficulties due to objective external 
factors and sometimes a pessimism as to purposefulness 
of developing the methods based on models (in the times 
of proliferation of computer simulation) the Centre in 
Iława today ranks as one of the best worldwide, equipped 
with the most suitable training models and water areas 
providing the widest training possibilities. Also, an unique 
approach to realization training tasks has been elaborated. 
As aresult the use of vacancies in training courses reaches 
100%, and representatives of many countries and organi- 
zations are turning to the Centre for advise or directly for 
realization of a similar undertaking on their territory. So 
far, as many as 2600 ship masters, first officers and pilots 
from almost 40 countries and all continents have taken 

part in various training courses in Iława. 
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